Saddle scars: Existence and applications 
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Abstract 

A quantum scar is a wave function which displays an high inten- 
sity in the region of a classical unstable periòdic orbit. Saddle scar are 
states related to the unstable harmònic motions along the stable man- 
ifold of a saddle point of the potential. Using a semiclassical method 
it is shown that, independently of the overall structure of the poten- 
tial, the local dynamics of the saddle point is sufhcient to insure the 
general existence of this type of scars and their factorized structure 
is obtained. Potentially useful situations are identified, where these 
states appear (directly or in disguise) and might be used for quantum 
control purposes. 

1 Introduction 



Until the early eighties it was widely believed that, for systems with an er- 
godic classical motion, the squared eigenfunctions must coincide, in the semi- 
classical limit, with the projection of the microcanonical phase space measure. 
This idea found sòlid ground on the mathematical results of Shnirelman[T|], 
Zelditch0 and Colin de Verdière0. On close scrutiny however, what these 
results state is that, for a quantum system that is classically ergodic, there is 
an eigenvalue sequence of density one such that the corresponding quantum 
densities ^(x)) 2 converge weakly to the Liouville measure. Therefore the 
observation of states that do not fit these expectations does not contradict 
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the mathematical results. For one thing the convergence may be very slow 
and, on the other hand, nothing forbids the existence of other subsequences 
converging to measures different from the Liouville measure. 

In fact wave functions were found which are concentrated near the classi- 
cal unstable periòdic orbits. When this happens one says that the quantum 
state is scarred by the unstable periòdic orbit or that one has a quantum scar. 
Such states have been observed at first in numerical simulationsQ] || || 
and, more recently, experimental evidence was found [0] on a semiconductor 
quantum- well tunneling experiment. 

The first theory of scars was proposed by Heller|5[], other theoretical for- 



mulations followed, developed by several authors|| [Kj []TTj] . Heller's 
theory studies the overlap integral 

C(t) = (*(t,x)\*(0,x)) (1) 

for a propagating wave packet which at time zero has a Gaussian shape 
and initial conditions (p ,x ) corresponding to an unstable periòdic orbit. 
Expanding \l/(0,x) in energy eigenstates 

^(0,x) = ^c„^ n (x) (2) 

n 

one sees that the Fourier transform S(E) of the overlap C(t) is the spectral 
density weighted by the probabilities |c n | . 

S(E) = Y J \c n \ 2 à{E-E n ) (3) 

n 

Now, if the period r of the classical periòdic orbit and the largest positive 
Lyapunov exponent A are such that e~ rA / 2 is not very small, the overlap C(t) 
will display peaks at times nr. As the wave packet spreads, the amplitude 
of the peaks decreases after each orbit traversal at the rate e _rA//2 . The 
Fourier transform of C(t) will therefore have peaks of width À with spacing 
uj = — . Referring to (|3|) one concludes that only the eigenstates that lie 
under the peaks contribute to the expansion of the wave packet. Since the 
wave packet has an enhanced intensity along the region of the period orbit, 
this is expected to carry over to the contributing energy eigenstates. The 
stronger the overlap resurgences are, the stronger the effect is expected to 
be. Therefore the intensity of the effect varies like 1/rÀ . 
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The above qualitative derivation|| of the scar effect is flawed if the prod- 
uct Xd(E) (where d(E) is the mean level density) is very large. Then the 
number of contributing eigenstates is very large and no individual eigenstate 
is required to show a significant intensity enhancement near the periòdic 
orbit. Also the argument assumes the low period unstable orbits to be iso- 
lated. If there are several nearby orbits of different periods the argument 
breaks down. However, if it happens that many periòdic orbits of the same 
period are present in the same configuration space region, the effect may even 
be enhanced. This is the situation for the periòdic motions in the neighbor- 
hood of an unstable critical point (a saddle) of the potential V(x) . Near 
the critical point there are unstable harmònic periòdic motions along the 
stable manifold of the critical point. As long as anharmonic corrections are 
unimportant, all the orbits will have the same period independently of their 
amplitude. The scars associated to these unstable periòdic orbits are called 
saddle scars in this paper. 

Whenever a dynamical systems has a phase-space region with sensitive 
dependence to initial conditions, the periòdic orbits in that region are unsta- 
ble and, even when they are dense, they are a zero measure set in the smooth 
(Liouville) measure over the energy surface. Therefore, unstable classical or- 
bits are in practice never observed, because all typical motions are aperiodic 
and uniformly cover the support of the Liouville measure. The phenomenon 
of quantum scars may therefore have far-reaching implications for the appli- 
cations of quantum systems. Whenever an unstable periòdic orbit scars a 
quantum eigenstate, the system may easily be made to behave like the un- 
stable orbit by resonant excitation to the corresponding energy level. In this 
sense, scars are a gift of Nature, for they allow the exploration of dynamical 
configurations that in classical mechanics are washed away by ergodicity. 

Saddle points are the typical critical points of genèric (Morse) functions. 
Therefore, once their existence is established, saddle scars are expected to 
be quite abundant. In the remainder of the paper a semiclassical method is 
used to establish that, independently of the overall structure of the potential, 
the local dynamics of the saddle point is suflicient to insure the existence of 
this type of scars. Then, in the closing section, potentially useful situations 
are identified where these states appear (directly or in disguise) and may be 
used for quantum control purposes. 
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2 Semiclassical estimates 



In the neighborhood of a saddle point, there is a choice of coordinates such 
that, up to higher order terms, the potential is 

V(x) = £ a t x\ + ■■■ (4) 

i 

For two dimensions o\ > and Oi < 0. In this case one obtains the following 
result: 

# There are scar states concentrated along the stable manifold 
of the saddle point and, on the neighborhood of the stable manifold, 

l^scarl 2 oc cos Ç^ x tj \Í>n{xi)\ 2 (5) 

v/j n (xi) being close to an harmònic oscillator wave-function and W n 
a function of the monodromy matrix in the transverse direction. 

(Explicit expressions for W n under different approximations are given below) 
The result is obtained following Bogomolny's semiclassical construction of 
wave functions[|J. From the series expansion of the energy Green's function 
in terms of eigenfunctions it follows that the averaged squared wave function 
is proportional to the imaginary part of the Green's function 

(|^ (x)| 2 ) oc (lmG{x,x,E )) (6) 

the average being taken over a small energy interval AE around E , which 
corresponds, in the semiclassical approximation, to restrict the contributions 
to orbits with times of motion of order < . Likewise, if an average is 
taken over small intervals of the variable x, the dominant contributions come 
from classical trajectòries for which the change of momentum on the closed 
orbit is small. I will discuss later the role of these two averages. 

The next step is to use the semiclassical approximation for the Green's 
function 

_ / 1 \ (d+l)/2 

G(x , x, E) = G(x , x, E) + í - J GqscOo, x, E) (7) 



Gosc{.x ,x 1 E) = 




(8) 
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For the contributions in the neighborhood of a classical periòdic orbit it is 
convenient to choose one of the coordinates along the orbit (xi) and the 
others (xj ; i = 2, ■ • ■ , n) along the transverse directions[T2|] ||. On the 
neighborhood of an unstable periòdic orbit along the stable manifold of the 
saddle point the action is expanded up to quadratic terms in the transversal 
coordinates and one has 



G OS c(x,x,E) 



1 



i (27TÍ 



,V2 



DV^Xi) \ i í 1 - 1 ,7T 

L, • ex P \ * \ S P + o 2^ WijixijXiXj | - i-up 



P 



i,3=2 



(9) 

with D(xi) and Wij(xi) functions of the monodromy matrix in the transverse 
coordinates 



with 



x±(n) \_í ^n(^i) ™>i2(xi) 

P±(n) J \ m 2 i(x 1 ) m 2 2(x 1 ) 

D(xi) = det (mï<f) 

(m 22 - l)mi 2 1 




(™T 2 ) 



(10) 



(11) 



W(xi) = m 12 mn 

and n the period of the periòdic orbit along x\. These expressions hold for 
any number of transverse directions. I now specialize to a two dimensional 
saddle point. Non-trivial orbits along the stable manifold are harmònic mo- 
tions of period T\ = 
Therefore defining 



2tïJ T^rj- independent of the amplitude of the oscillation. 



D 



^2m 2 a 2 



sinh ( 2vrJ^^ 

rri2<T\ 



2 cosh 2tt 



w 



m 2 cri 



- 2 



■ sinh 2tt 



(12) 



(13) 



v / 2m2 0-2 

D and W do not depend on xi and the dependence on the transverse co- 
ordinate factors out for each term of the sum in Eq.(^). However, for each 
primitive trajectory, one also has to sum over múltiple passings obtaining the 
following sum 



(2m 2 a 2 ) 4 sinh 2 (n8) exp 




(14) 
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where 9 = 2n./ :ni£ ^ . There are two situations where simple closed form 
results may be obtained: 

# When exp (6>) ^> 1 , the a^-dependence factors out from the sum and 
the result is 

(\^e(x)\ 2 ) oc (imjexp G(x u E)}} (15) 

G(x\, E) being the harmònic oscillator Green's function. From 

g( Xi , e)=y: \* n (x)\ 2 [p - m5 ( E - ( ie ) 

it follows that the principal part drops out under averages over small energy 
intervals and the result (0) follows with W n = W for n = 1, 2, • • • . The lowest 
lying state however corresponds to the orbit of the unstable fixed point and 
the period of this orbit is no longer r%. The value of Wq in this case may 
be calculated by considering an orbit from the coordinate (0, £2) to the fixed 
point (0, 0) and back along the unstable manifold, which is equivalent to take 
the limit a% — > in Eq. (|Ï3|) . Then 

W = 2 v / 2^ 2 " (17) 

# If exp (9) is not large then the 22-dependence does not factor out in 
the sum flï~4]) , Notice however that the sum in (11) is not a sum over múltiple 
passings of the same orbit, but a sum over different orbits because, for x 2 7^ 
, the initial and final momentum are different, and the difference grows with 
n. Therefore if, in addition to the average over small energy intervals, one 
also averages over small coordinate intervals then, the contribution of the 
primitive orbit dominates the leading semiclassical approximation. The same 
factorization of the X2-dependence, as before, is obtained. However, the ip{x{) 
function in this case is not exactly an harmònic oscillator wave functions, but 
a function corresponding to a sum restricted to the primitive orbits. 

3 Applications 

The canonical form of the potential near a saddle point establishes a local 
separation of variables which, of course, does not hold far away from the sad- 
dle point. However, what the semiclassical estimate of the previous section 
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shows is that the local dynamics of the saddle point is sufíicient to insure 
the local existence of a factorized quantum state fl5|). If the separation of 
variables extends over a sufficiently large range tlien, the transverse shape of 
the wave function may be approximated by the solution of a one dimensional 
problem. Consider, for example a one-dimensional Hamiltonian 

H= -éi& + 2scoa< - x) (18) 

with 2tc— periòdic boundary conditions. The eigenstates are Mathieu func- 
tions and there is a state of energy slightly above 2g with the squared am- 
plitude |F| 2 as shown in Fig.l (for g = 50 , ^ = 1 , E = 101.189). In 
the figure the state is also compared with the local approximation following 
from (U) and (|T7|). The width of the state depends on the factor ^/2mg and 
the energy is approximately 2g plus the kinetic localization energy. For an 
higher dimensional problem one must also add the quantized energy of the 
harmònic oscillation along the stable manifold. These considerations provide 
a simple rule for the approximate energies at which saddle scars are expected 
to be found. 

Among the physical situations where saddle scars might appear, an inter- 
esting example is probably the quantum collision of systems containing both 
attractive and repulsive interactions (chemical ions, nuclei, etc)|]T3[. When 
a system contains several positive and negatively charged partides, there 
are classical configurations of close proximity of the partides which are of 
low energy because the repulsion between the like-charged partides is com- 
pensated by the attraction of unlike-charged partides. These configurations 
however are highly unstable and the chance to observe (or stabilize) them in 
classical mechanics is nil. They are zero measure configurations in the en- 
ergy surface. For smooth potentials these unstable configurations would be 
saddle points of the potential, hence they are expected to give rise to saddle 
scars. These states would correspond to well defined energy levels and might 
be prepared by resonant excitation. That is, quantum control through scars 
makes accessible some states that, classically, are essentially unobservable. 

Another interesting aspect of saddle scars is their generality, because 
saddle points are the typical critical points of genèric functions. There are 
also other features of the classical phase space which the wave functions 
imitate and that, in some limit or by change of coordinates, may be related 
to the saddle scars. This concerns in particular the regularization of singular 
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potentials. An example is the collision states found for a three-dimensional 
periòdic Coulomb problem|13]. In this case the quantum collision states 



correspond to wave functions concentrated along a phase space feature which 
is not an actual orbit, but the separatrix of two classes of unstable orbits. 

In Jacobi coordinates (r = X\ — x 2 , r] = x 3 — | {x\ +2:2)) the potential 
between three unlike-charged partides is 



r 



-rj 



+ T) 



(19) 



The dynamics of binary col·lisions in the three-body problem may be regu- 
larized, but the case of interest here is a triple collision which, except for 
exceptional cases[]nj, is not regularizable. The potential, however, may be 
regularized by addition of a small quantity to the definit ion of the distances 



H e = (pí + PÍ + PÍ + e' 



In the neighborhood of the triple collision point 
potential is 



77 



(20) 

, the regularized 



V f = 



1 1- 

\v\ + ■ 



(21) 



1 I 1 

— r H 

e 4e 3 e 3 

and this is a saddle point with three stable and three unstable directions. 
According to the discussion above one would expect the quantum collision 
states to have scarred wave functions concentrated along the stable manifolds 
and with a small dispersion in the transverse (unstable manifold) directions. 
This is the situation that is indeed found in the numerical computations[|T3| 
of the three-dimensional periòdic Coulomb problem. It shows that the ef- 
fect seems to survive the e — >• limit. Triple col·lisions in a 3-dimensional 
3-body problem lie on an analytic 10-dimensional submanifold of the 12- 
dimensional dynamical manifold. Hence it is a zero measure effect, essen- 
tially non-observable in classical systems. It is interesting that, through the 
scar effect, they do correspond to well-defmed energy levels, accessible by 
resonant excitation. 



4 Figure caption 

Fig.l - One dimensional density for a wave function concentrated around an 
unstable point and the semiclassical approximation (+). 
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